GW approximations and vertex corrections on the Keldysh time-loop contour: 
application for model systems at equilibrium 
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We study the effects of self-consistency and vertex corrections on different GW^-based approxima- 
tions for model systems of interacting electrons. For dealing with the most general case, we use the 
Keldysh time-loop contour formalism to evaluate the single-particle Green's functions. We provide 
the formal extension of Hedin's GW equations for the Green's function in the Keldysh formalism. 
We show an application of our formalism to the plasmon model of a core electron within the plasmon- 
pole approximation. We study in detail the effects of the diagrammatic perturbation expansion of 
the core-electron/plasmon coupling on the spectral functions in the so-called S-model. The S-model 
provides an exact solution at equilibrium for comparison with the diagrammatic expansion of the 
interaction. We show that self-consistency is essential in GW-based calculations to obtain the full 
spectral information. The second-order exchange diagram (i.e. a vertex correction) is also crucial 
to obtain the good spectral description of the plasmon satellites. We corroborate these results by 
considering conventional equilibrium GW^-based calculations for the pure jellium model. We find 
that with no second-order vertex correction, one cannot obtain the full set of plasmon side-band 
resonances. We also discuss in detail the formal expression of the Dyson equations obtained for the 
time-ordered Green's function at zero and finite temperature from the Keldysh formalism and from 
conventional equilibrium many-body perturbation theory. 
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I. INTRODUCTION 

Equilibrium, zero- and finite-temperature Green's 
functions techniques based on many-body perturbation 
theory (MBPT) are widely used in electronic-structure 
and total energy calculations^. Hedin's formulation^ 3 - 
for the electronic Green's function closes the many-body 
hierarchy by expanding the electron self-energy of the 
one-particle Green's function in terms of the screened 
Coulomb interaction in the presence of vertex corrections. 

Without these vertex corrections, one obtains the con- 
ventional GW equations 3 - - — . The GW method is an ap- 
proximate treatment of the propagation of electrons: it 
can be seen as if electrons interact with themselves via a 
Coulomb interaction that is screened by virtual electron- 
hole pairs. In bulk semiconductors, the GW approxi- 
mation is known to lead to surprisingly accurate band 
gaps^i&iiS, while for finite-size systems and molecules the 
method provides qualitatively correct values of ionization 
energies and electron affinities^. It also provides a con- 
venient starting point for many useful approximations 
and applications to photoemission spectroscopy^ and op- 
tical absorption in metals or semiconductors as well as in 
finite size molecular systems 6 -^ - — . Most practical GW 
calculations today are performed in a perturbative man- 
ner using equilibrium MBPT. 

However, if we want to consider a system driven out 
of equilibrium by an external "force" , such as, for ex- 
ample, a molecular wire coupled to electrodes sustaining 
an electronic current flow, or any system driven by an 



external electromagnetic field (time-dependent or not), 
we need to extend the equations for the dynamics of the 
quantum many-body interacting system (Hedin's equa- 
tions or their simplified GW form) to non-equilibrium 
conditions. 

For this, the non-equilibrium Green's function (NEGF) 
technique^— has been widely used to calculate 
electronic transport properties of mesoscopio« and 
nanoscalei^— systems, plasmas, quantum transport 
in semiconductors^ and high-energy processes in nu- 
clear physics 2 ^. Also known as the closed time-path 
formalis m 26 ' 27 , the NEGF formalism depends on an "ar- 
tificial" time parameter that runs on a mathematically 
convenient time-loop contour (plus eventually an imagi- 
nary time for taking into account the initial correlation 
and statistical boundary conditions). It is a formal pro- 
cedure that only has a direct physical meaning when one 
projects back the time parameters of the time-loop con- 
tour onto real times. It was introduced because it allows 
one to obtain self-consistent Dyson-like equations for the 
Keldysh Green's function using Schwinger's functional 
derivative technique. Transforming the Dyson equation 
to real time by varying the Keldysh time parameter over 
the time-loop contour results in a set of self-consistent 
equations for the different non-equilibrium Green's func- 
tions (advanced/retarded or lesser /greater). The NEGF 
technique is general and can treat non-equilibrium as well 
as equilibrium conditions, and the zero and finite temper- 
ature limits, within a single framework. 

The NEGF technique has been applied to the study 
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of different levels of self-consistency in the GW a ppr oach 
for atoms, molecules and semiconductors in Refs. [28l - l34j . 
However, these works did not include the effects of simul- 
taneous self-consistency and vertex corrections. Other 
levels of approximation for electron-electron interactions 
have also been considered in finite-size nanoclusters by 
using the Kadanoff-Baym flavour of NEGF^>2&. 

In this paper, we want to study these effects (self- 
consistency and vertex corrections) and use the most gen- 
eral formalism to deal with the full equivalent to Hedin's 
GW equations. We believe that the Keldysh formalism, 
even applied to equilibrium conditions, can be more use- 
ful than the conventional approaches since it is by nature 
a more general approach. 

We extend Hedin's equations to the Keldysh time-loop 
contour, and derive the equations for the one-particle 
Green's function G, self-energy E, screened Coulomb in- 
teraction W and for the (3-point) vertex functions T. 
Note that a non-equilibrium approach to Hedin's GW 
equations has been provided in Ref. [37} where an al- 
ternative distinct approach based on the Liouvillian su- 
peroperator formalism is used. However, working in a 
Louivillian vector space is much less convenient and much 
more computationaly demanding for practical applica- 
tions than working within an Hilbert space as in the for- 
malism we develop below. 

We then apply our formalism to the calculation of 
the spectral function of a particular model of an ho- 
mogeneous electron gas: the plasmon model for a core 
electron 3 - ' 38 ' 39 . We choose this model as it can be solved 
exactly at equilibrium, and thus we are able to com- 
pare the different approximations introduced in the cal- 
culations (self-consistency versus one-shot calculations, 
and/or vertex corrections) and check their validity for 
different limiting cases (the high and low electronic den- 
sity regimes). We also compare the outcome of these 
calculations with conventional GW calculations for the 
jcllium model. We examine if the effects on the spectral 
functions rendered by self-consistency iterations and the 
inclusion of vertex corrections which we find for the plas- 
mon model with a core electron also hold for the jellium 
model. 

To our knowledge, the only available exact results are 
for equilibrium conditions, and thus we benchmark our 
formalism against exact results at equilibrium before ex- 
tending the discussion to non-equilibrium conditions. 

The paper is organized as follows. In Section|Hl we re- 
call the expressions of Hedin's GW equations and briefly 
review the performance of conventional GW calculations. 
The extension of Hedin's equations to the Keldysh time- 
loop contour is provided in Section !!!!] We also show that 
we recover the conventional non-equilibrium GW formal- 
ism developed and used by other s 28 ' 31 " — when ignoring 
the vertex corrections in Appendix [C]) . The lowest-order 
expansion, in terms of the interaction for the screened 
Coulomb interaction W and for the vertex functions T, 
is given in Appendix [Cj In Appendix [D] we also provide 
a rigorous mathematical proof of the difference between 



equilibrium time-ordered Green's functions in the zero 
and finite temperature limits that were discussed less rig- 
orously in Chapter IV. 17 of Ref. 0]. 

In Section IIV1 we apply our formalism to the calcula- 
tion of the spectral function of a model system and core 
electron coupled to a plasmon mode 3 - ' 38 ! 39 . The exact so- 
lution of this model at equilibrium permits us to exam- 
ine the effects of self-consistency and vertex corrections 
on the spectral density. We also examine these effects 
for another model of an electron gas, the jellium model, 
by using conventional GW calculations fSection IIVEI) . 
We show a general trend: second-order diagrams for the 
interactions (i.e. vertex corrections) are necessary to ob- 
tain the full series of plasmon side-band peaks. Finally, 
we conclude our work in Section [Vj 



II. HEDIN'S GW EQUATIONS 

Hedin's GW equations 2 -* 3 - were originally derived for 
the time-ordered single-particle Green's function G at 
equilibrium, defined by 



G(12) = -i(T*(l)* t (2)) 
They are expressed as follows^ 3 -: 



(1) 



G(12) = G (12) + / d(34) G (13) £(34) G(42), (2a) 



(2b) 



£(12) = i J d(34) G(13) T(32;4) W(41), 



P(12) = -i J d(34) G(13) G(41) T(34;2), 



r(12;3) = 5(12)5(13)- 



W(12) = v(l2) + I d(34) u(13) P(34) W(42), (2c) 

(2d) 
(2e) 

d(4567) G(46) G(75) T(67;3), 

with the usual notation for the space-time coordinates: 
any integer i represents a point in space-time i = Xi = 
(rj,tj); and for the electron single-particle Green's func- 
tion G, the corresponding self-energy £, the screened 
Coulomb interaction W, the irreducible polarizability P 
(sometimes also called polarization), and the vertex func- 
tion r. 

Up to now, most practical GW calculations are per- 
formed not fully self-consistently, using a single itera- 
tion of the GW equations, called one-shot GW or GqWq. 
When a single iteration is performed, the initial approxi- 
mation must be good, so typically Go is constructed from 
the orbitals of any current ground-state method which 
correctly predicts the basic physics of the system. 

The application of Go Wo corrections to spectral prop- 
erties and band gaps as calculated in the local density ap- 
proximation (LDA) and generalized gradient approxima- 
tion (GGA) in density functional theory (DFT) is a long- 
standing success story^S, at least for many s-p bonded 
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systems. However, total energies calculated from the 
Galitskii-Migdal formula at the GqWq level are generally 
worse than those given by other ground-state methods. 
The initial close agreements with measured band gaps 
were later shown to be partly due to technical approxima- 
tions used along the way. Several studies have shown that 
LDA+GoWo systematically underestimates band gaps 
when solved in state-of-the-art all-electron schemes with 
full explicit treatment of frequency integrals (gaps are 
underestimated by about 1-10% for s-p bonded systems 
and about 20-50% for systems with d electrons, like rare- 
earth oxides, sulfides and nitrides) 41 ' 42 . The remaining 
discrepancy has prompted the search for more accurate 
but still tractable methods. 

In general, attempts at fully self-consistent GW have 
shown that spectral properties worsen as compared with 
Go Wo, while total energies improve. GW band gaps 
were first shown to be larger than expected in a quasi 
one-dimensional Si wire model by de Groot et al^. 
Von Barth and Holm showed that in jellium for a self- 
consistent update only of the Green's function (a GWo 
approach) a displacement of weight from quasiparticle 
peaks into the incoherent background occurs 4 ^. Also, 
the occupied bandwidth broadens rather than narrows, 
as expected from experiments on simple metals. At the 
same time, E. L. Shirley showed that the bandwidth of 
jellium broadens further with full self-consistency, while 
the bandwidth narrows again once vertex corrections are 
taken into account 4 ^. The effects of non-locality in vertex 
corrections were also addressed in Ref.[46j. 

Later studies have shown that GW total energies for 
jellium are very accurate 4 ^—, as expected from a con- 
serving approximation in the Baym-Kadanoff sensed. 
This holds true even for low-dimensional atomic and 
molecular systems, and ionization potentials as calcu- 
lated by the extended Koopman's theorem also tend to 
be accurat o 29 ' 34 ' 51 . Full GW calculations were performed 
by Kutepov et al. for simple metals and semiconductors, 
They showed inter alia that the calculated equilibrium 
lattice parameters were all very close to the experimen- 
tal ones^ 2 -. 

In view of improving the starting point, quasi-particle 
self-consistent GW has emerged as a good compromise 
between self-consistency and a practical path to good 
spectral properties^ 3 -. It has been shown that vertex cor- 
rections further improve the correspondence between the- 
ory and experiment^, but consistently accurate results 
still remain elusive for systems with localized states, de- 
fects and band offsets^— . 

In particular relevance to the present paper, no exis- 
tant implementation of GW seems to describe the full 
spectrum of plasmon satellites in metals. 



III. EXTENSION OF HEDIN'S GW EQUATIONS 
TO THE KELDYSH TIME-LOOP CONTOUR 

We now consider the generalization of the single- 
particle Green's function on the time-loop contour (the 
so-called Keldysh contour Ck with two branches, branch 
(+) for forward time evolution and branch (— ) for back- 
ward time evolution) : 

G(12) = -i(T CA .*(l)* t (2)) . (3) 

For the moment, we do not specify the nature of the 
"external force" that drives the system out of equilib- 
rium. We consider the generalized Green's function on 
the Keldysh time-loop contour and hence end up with 
four different Keldysh components for the Green's func- 
tions: G ++ ,G H ,G~ + ,G , defined according to the 
way the two real-time arguments (ti,^) are positioned 
on the time- loop contour Ck- The initial correlations (i.e 
the initial boundary conditions) are assumed to be dealt 
with in an appropriate wa y 15 ' 17 ' 22 . 

To derive the NE-GW equations, we proceed as fol- 
lows: in each integral Jd(l), the time is integrated 
over the time-loop contour Ck- J c dr% and then de- 
composed onto the two real-time branches: j c dr% = 
J, + , dtf + J*,_n dt^ = J di+ — J dt±. We then calculate 
the different components X mV2 (with r/i2 = ±) for the 
Green's function, self-energy, screened Coulomb interac- 
tion W, polarizability P, and vertex function T. Where 
possible, we re-express these in a more convenient way 
by using the relations between the different Green's func- 
tions and self-energies on the time-loop contour (see Ap- 
pendix |2| ■ 

There are actually three kinds of equation in Hedin's 
GW equations Eq. (2a-e). First, there is a set of Dyson- 
like equations for the electron Green's function G and for 
the boson Green's function W, i.e. the screened Coulomb 
interaction. In these two equations, the vertex func- 
tion r does not appear explicitly. Next there is another 
set of equations for the electron self-energy E and for 
the polarizability (the boson self-energy) P. In these 
equations, the vertex function appears explicitly. Fi- 
nally there is the equation for the vertex function it- 
self, r. The vertex function can be expanded as a series 
r(12;3) = Yin ^(n)(12; 3) where the index n represents 
the number of times the screened Coulomb interaction 
W appears explicitly in the series expansion. Each oc- 
currence of the screened Coulomb interaction W in the 
vertex function T is generated by the functional deriva- 
tive SE/6G. 

Finally, one should note that the equilibrium proper- 
ties of the system are, in principle, recovered from the 
extension of Hedin's GW equations to the Keldysh time- 
loop contour when the external driving force is omitted 
and the whole system is at thermodynamical equilibrium. 
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A. The electron Green's function and the 
self-energy 

Following the prescriptions given above, we calculate 
the components G ++ , G H and G h from the exten- 
sion of Eq. ((2j on the time-loop contour, and we find the 
Dyson-like equation for G r,a : 



G r ' a (12) = G r ' a (l2) + J d(34) G£' a (13)£ r ' a (34)G r ' Q (42), 

(4) 

which has the same functional form as in Eq. @. 

We also obtain the following quantum kinetic equation 
(QKE)forG^: 

G^(12) = J d(3456) 

[5(U) + G r (13)S7(34)] Gf (45) [5(52) + £ a (56)G Q (62)] 
d(34) G r (13) £^(34) G a (42). 

(5) 

B. The screened Coulomb potential 

By looking at Eq.(2c), one can see that W has the 
same functional form as the electron Green's function G. 
The screened Coulomb interaction W is a bosonic Green's 
function with an associated bosonic self-energy, the po- 
larizability P. With the formal equivalence (G, E) -H> 
( W, P) , one can expect to obtain a Dyson- like equation 
for the advanced and retarded screened Coulomb inter- 
action and a quantum kinetic equation for as equiv- 
alently obtained for the electron Green's function. 

This is indeed what we find: W r ' a follows the usual 
Dyson-like equation as 



W r ' a (12) = v(12) + I d(34) W (13) P r < a (34) W r < a (42), 

(6) 



or m a more compact notation 



W r ' a = v + vP r ' a W r ' a = v + W r ' a P r ' a v 
= v[l - P^v}- 1 = [1 - vP'-^v 



(7) 



where any product XY implies a space-time integration 
[Xy](12)= fd(3)X(13)Y(32). 

Since the bare Coulomb potential v(12) is instanta- 
neous, it corresponds to an interaction local in time and 
therefore its extension to the Keldysh contour has no 
it or v ^ components. Hence, we obtain the following 
quantum kinetic equations for : 

W^{12) = [ d(34) W r (13) P^(34) W a (42). (8) 



C. The vertex function T(12; 3) on the contour Ck 

The derivation of T(12; 3) on Ck does not create any 
formal difficulties. However since T(12; 3) is a three-point 
function, it is not possible to recover a Dyson-like or a 
quantum- kinetic- like equation for T. 

For any Keldysh components of the vertex function 
r V3V2V4 (32; 4) , we can formally write the different com- 
ponents of the self-energy on the Keldysh contour as fol- 
lows: 



£"i" 2 (12)=i^ mm J d(34) 



G mr,3 (l3) T I ' 3 ' ?2 ' )4 (32;4) W mvi (41), 



(9) 



and likewise for the polarizability 



p^(12) = -i^ mm / d(34) 

»)3')4 

G mr,3 {13) C 4 ' )1 (41) T' )3 ' )4 '' 2 (34;2). 

(10) 

Now we need to close the above equations, i.e. to find 
an equation for the different components Y Vin2m (12; 3) 
of the vertex function. By considering the equivalent of 
Eq.(2e) on the Keldysh contour, we obtain 

F&»»»»(i2;3) = ^ i " 2 (i2)5" 1 ' ?3 (i3)+ mvsmm 

)) 4 ...»)7 



d(4567) 



SGw* (45) 



G mve (4Q)G V7V5 (75)^^3 ( 67; 3 ) 

(11) 



In Appendix O we consider the series expansion of 
the vertex function T(12;3) = J2 n (12; 3) where the 
index n represents the number of times the screened 
Coulomb interaction W appears explicitly in the series 
expansion, and we provide explicit results for the elec- 
tron self-energy £ and polarizability P for the lowest or- 
der terms r (0) (12; 3) and r (1) (12; 3). 



IV. APPLICATION TO MODELS RELATED TO 
THE HOMOGENEOUS ELECTRON GAS 

Now we want to test our extended formalism of Hedin's 
GW equation onto the Keldysh time-loop contour and 
the corresponding series expansion of the vertex func- 
tions. The importance of self-consistency and vertex cor- 
rections was discussed in Section [TTJ Self-consistency and 
vertex corrections apply in both equilibrium and non- 
equilibrium systems and therefore are more conveniently 
addressed in as simple a model system as possible. 

Calculations could be performed for several model sys- 
tems, but would not lead to any pertinent conclusions 
if they could not be compared to exact results. To our 
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knowledge, exact results for interacting electron systems 
are few and not as widespread as numerical (highly accu- 
rate) calculations even for models of interacting electron 
systems. One of the available exactly-solvable models has 
been used in the context of x-ray spectroscopy of metals, 
and leads to tractable analytical expressions for the elec- 
tron Green's function: the plasmon model for the core 
electron^. 

In the next section, we consider this exactly-solvable 
model and compare the exact results with those obtained 
from our GW formalism, at zero and finite temperatures 
and with or without lowest-order vertex corrections. We 
note here that the exact solution is obtained for a model 
of an homogeneous electron gas at equilibrium. Dealing 
with an interacting system at equilibrium does not cause 
any problem within our formalism, since the equilibrium 
condition is just a special case of our more general for- 
malism for non-equilibrium conditions (See appendix [D] 
for a full discussion about the equilibrium limit of the 
Keldysh formalism at zero and finite temperatures). 



persion: 
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H ) 2 3 ( W ° 



0\2 



1/2 



(14) 



with uj° p = Lu p /e F = 4(^a)V2, a = {±Y'\ and r s 
defines the electron density n = (-^l) - " 1 



B. The S-model 



A particularly simple model of a core electron, known 
as the S-model 3 -^, is obtained by further replacing w" 1 



by a step function u> 



-l 



1 9(q c — q), where the cut-off 



parameter q c is determined by: 



dq 



2 d 1- 



(15) 



From this definition of q c , it follows that the energy shift 
parameter 



Effective Hamiltonian for the plasmon model of 
a core electron 



q ■ 



(16) 



The properties of an homogeneous 3D electron gas can 
be well-described within the plasmon model. The plas- 
mon model is defined from Hedin's equations Eqs. (2a- 
e) together with the so-called plasmon-pole parametriza- 
tion. In reciprocal space, the screened Coulomb poten- 
tial can be written as W(w, q) — v q e _1 (u;, q), where v q is 
the Fourier component q of the Coulomb potential. The 
dielectric function e _1 (w,g) is then obtained from the 
plasmon-pole approximation e -1 ^,^) = 1 + u> 2 i /(lu 2 — 
uiq), where uj p is the bulk plasmon energy, related to the 
electron density n as usual, uj 2 = (Airne 2 /m), and the 
plasmon dispersion u q remains to be defined. 

Within this model, the dynamic part of the Coulomb 
potential W(u),q) — v q can be re-expressed as 



v 2 = v q (e 1 (u,q) - l) = 



2uj„ 



2uj n 



= l q B(u,q), 

which involves a coupling constant j q and the bosonic 
propagator B(u,q) of the plasmon modes. 

Following Refs. [3.38.39] we consider the following 
Hamiltonian for the plasmon model of a core electron 



e c c^c + y]u) q blb q 



J2l^c(b q + bl q ). (13) 



For this model of the core-electron case there exists a 
precise and well-defined relation between the solution de- 
fined by a plasmon model for an electron gas and the 
solution defined by the corresponding effective Hamil- 
tonian H e ff^. Finally we consider the q — > limit of 
static random-phase approximation 3 for the plasmon dis- 



is the same as for the corresponding plasmon model. 

The solution of the S-model can be mapped onto a 
simpler Hamiltonian, giving rise to the same spectral in- 
formation 

H cS = eje + Lo p tfb + 7oc f c(6 + 6 1 ), (17) 

with 7g = Dujp. 

An analytical expression for the relaxation energy D is 
found from the chosen dispersion relation of the plasmon 
frequency uj q . 

We then find that the corresponding relaxation energy 
is given by 



^oj„ 2 J (2irY Vq ^ a 2 V 



2^K + f)V2- 

(18) 

This result is very similar to the relaxation energy 
found by Minnhagen 3 ^ when one replaces the prefac- 
tor 4r in the dispersion relation w q by I and when one 



uses the trigonometric relations sin(^) = 



and 



cos[taa _1 («)] = 1/V(1 + u 2 ). 

The other advantage of dealing with the S-model 
is that it has an exact solutio n 38 i 39 i 58 which can be 
compared with approximate calculations performed with 
Hedin's GW equation for different levels of expansion of 
the self-energy and/or vertex function. The exact solu- 
tion of the S-model at zero temperature provides us with 
an analytical expression for the retarded Green's func- 
tion, given by 



y 2 7 



2d 



1 



n=0 



n\ uj — e c + nujp + \rj 



(19) 
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with 7 2 = (70/wp) 2 — D/tUp and the renormalized core 
level s c = e c + D = e c + j 2 ujp. The finite tempera- 
tures solution is obtained from the prescription given in 
Ref. 



C. Feynman diagrams for the self-energy 

The Hamiltonian for the S-model given by Eq. p7|) is 
effectively a single electron coupled to a single-boson- 
mode model similar to the model we studied for an 
electron-phonon coupled system in Refs. 59.60]. We can 
then use the NEGF code we have developed to study the 
electronic properties of the S-modcl for different levels of 
approximation for the corresponding self-energies. In the 
Feynman diagram language, these are given in Figure [T] 
and correspond to (a) non-self-consistent calculations for 
the self-energy E = GqW P7 where Go is the core-electron 
bare Green's function and W p is the plasmon propagator 
given in Eq. (TT2"1) ; (b) self-consistent calculations for the 
core electron Green's function E = GW P ; and to vertex 
corrections taken at the Tm level of approximation for 
(e) non-self-consistent calculations E = GT^^W V with 
G and TW taken at the GW P level and (f) fully self- 
consistent E = GT^Wp calculations. 

Our NEGF code, presented in Ref. 59] is versatile. It 
was originally developed to deal with an electron-phonon 
coupled system in contact with two electron reservoirs 
each at their own equilibrium. But the code can deal 
with any model Hamiltonian of electron-boson coupled 
systems. In the following we use this code and we con- 
sider the whole system at equilibrium, and at zero or fi- 
nite temperature. As explained above, the exact solution 
of the S-model exists only for the equilibrium condition. 

Additionally we use an extremely small coupling con- 
stant to the reservoirs in order to introduce a finite but 
very small broadening in the spectral features of the S- 
model Hamiltonian Eq. (jTTJ) in a simple way (77 has a tiny 
but finite numerical value). The details for the calcula- 
tions of the different NEGF, at equilibrium and out of 
equilibrium, are given in Ref. [59l| . 

In Ref. [59], we discussed the first and second- 
order diagrams for the electron-phonon interaction — 
topologically speaking, this will look similar to the GW- 
like self-energy diagrams we consider here (Fig.[TJ, how- 
ever there the boson line is the phonon propagator and 
not the screened Coulomb interaction W with which we 
are concerned here. Furthermore the parameters of the 
core electron-plasmon coupled system are given here by a 
single physical quantity: the electron density (see Table 



D. Results 

Within our model, all the characteristics of the plas- 
mon are determined by a single parameter: the electron 




FIG. 1: Different levels of approximation for the one-particle 
self-energy E = E' 1 ' + E' 2 ' within the plasmon model. First 
order diagrams: (a) E^ 1 ' = GqW p with Go being the bare 
core-electron Green's function; (b) E^ 1 ' = GW V for self- 
consistent calculations. Second order diagrams with vertex 
corrections (c) E 1 - 2 ' = GTf^Wp for non self-consistent calcu- 
lations; (d) E (2) = G rf^j W p for the full self-consistent calcu- 
lations (see Appendix IC2|) . 
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0.02984 


EF 


0.0737 
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0.2046 


0.4604 


UJp 


0.1549 


0.2165 


0.3333 


0.6124 




2.103 


1.881 


1.629 


1.330 


D 


0.34046 


0.31186 


0.27789 


0.23523 


70 


0.22966 


0.25985 


0.30435 


0.37953 


70 /Up 


1.48 


1.20 


0.91 


0.62 



TABLE I: Values (in atomic units) of the different relevant 
parameters, electron density n, Fermi energy ef, plasmon en- 
ergy uj p , electron-plasmon coupling constant 70 and relaxation 
energy D for different values of rs- 

density, or equivalently by the Wigner-Seitz radius r$. 
There is then only one other parameter left: the energy 
level e c of the core electron, which we take as being lo- 
cated one atomic unit of energy below the Fermi level Ef 
of the different systems we consider. 

Table Q] contains the values of the different relevant 
parameters for four different values of rs- The high- 
density limit (rs = 2) corresponds to a medium electron- 
plasmon coupling, while the low-density limit (rs = 5) 
corresponds to a very strong electron-plasmon coupling. 

Below, and in Figs. [2HSJ we show results for the spec- 
tral function A(uj) = -i(G r (w) - G a (u))/2 calculated 
at equilibrium for two values of rs (medium coupling 
rs = 2, strong coupling rs = 4) at zero temperature 
and at a finite temperature. We compare the exact re- 
sults Eq. (fT9"| for the spectral function with the results 
obtained from the diagrammatic expansion of the self- 
energy and the vertex function shown in Fig. [TJ 

1. Exact results 

The exact spectral function, calculated from the ex- 
pression for the Green's function given in Eq. (|19p . is 
shown as a solid black line in Figs. [2] and [3l A broad- 
ening equal to the broadening of our NEGF calculations 
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FIG. 2: (Color online) Zero-temperature equilibrium spectral 
functions A(u>) for the for the high-density limit with rs = 
2, corresponding to medium core electron-plasmon coupling 
"fo/ui p — 0.62. Top panel: Exact results and GW calculations 
with and without self-consistency E = GW P , GoW p . Bottom 
panel: Results for different levels of approximation for the 
self-energy E = G W P , GW P , G(F (0) +T^)W P and G(r (0) + 
Tf^)W p (see Fig. [TJ with fewer grid points (N w = 1579), 
giving an extra broadening. 




(co - eJ / co 

v F 7 p 



FIG. 3: (Color online) Zero-temperature equilibrium spec- 
tral functions A(uj) for the low-density limit with rs = 4, 
corresponding to very strong core electron-plasmon coupling 
7o/oj p = 1.20. Top panel: Exact results and GW calculations 
for the different self-energies E = GoW p , GW P . Bottom panel: 
Results for different levels of approximation for the self-energy 
E = G W P , GW P , G(T( ) + T$f)W P and G{T {0) + T$)W P 
(see Fig. [TJ with fewer grid points (JV W = 1579), giving an 
extra broadening, in comparison to the top panel. 



has been applied. Fig. [5] shows the zero-temperature re- 
sults for the high-density electron gas (rs — 2). The 
exact result provided by Eq. (fTI?)) (solid black line) gives 
a spectral function with a peak localized at the renormal- 
ized core level e c = e c + D, and plasmon side-band peaks 
at e c — nuj p (n > 1) corresponding to plasmon emission. 
The peaks are hence separated by the plasmon energy uj p . 
In terms of amplitude, the main peak is that at e c in the 
limit of weak to medium/strong electron-plasmon cou- 
pling, i.e. where jo/u> p < 1, and so for which 7o/w p < 1. 

Fig. [3] shows the zero-temperature results for stronger 
coupling, rs — 4 and 70/wp > 1- Now the renormal- 



ized core level ec has shifted to oj/u> p ~ 1.5, while the 
spectral weight is shifted towards lower energies and the 
main peak is now the plasmon side-band peak at around 

2. Diagrammatic expansion results 

The main differences between the exact result and the 
diagrammatic expansions of the self-energies and of the 
vertex functions (as represented in Fig.[TJ are as follows: 

First, let us discuss the results for the spectral func- 
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tions in the high-density limit (rs — 2) for which the 
electron-plasmon coupling is medium "fa/tOp = 0.62. 

The non-self-consistent GW calculations (i.e. S = 
GoW p , Fig. HJa), dotted black lines in Fig. [5J generate 
only two peaks, the renormalized core level with one plas- 
mon side-band peak, as expected. However the positions 
of those two peaks are incorrect. 

The self-consistent GW calculations (i.e. £ ~ GW P , 
Fig. [2(b), solid green lines in Figs. [5] and [3]) generate the 
correct series of plasmon side-band peaks. However the 
corresponding relaxation energy D is too small and the 
energy position of the first plasmon side-band peak is too 
low. It should be noticed however that the energy sepa- 
ration between the plasmon side-band peaks is correctly 
reproduced, i.e. equal to ui p . 

For the low-density limit (rs — 4) for which the 
electron-plasmon coupling is very strong jo/cu p = 1.20, 
the GW calculations poorly describe the exact spectral 
density. The self-consistent GW calculations generate 
the correct series of peaks but with a completely wrong 
weight distribution. This is unsurprising since the GW 
approach corresponds to a partial resummation of the di- 
agrams, and does not include all other relevant diagrams 
necessary to deal with the very strong regime. 

The lowest-order vertex corrections to the self-energy 
(Figs. [TJd) and [He), blue dashed lines and red trian- 
gles in Figs. [5] and [3]) introduce modifications of the peak 
positions. They generate a slightly better relaxation en- 
ergy D and a shift of the side-band peaks towards the 
renormalized electron core level (Figs. [2] and [3j bottom 
panels). Vertex corrections globally improve the spec- 
tral information towards better overall agreement with 
the exact results. However, the lowest-order vertex cor- 
rection expansion F( ) + Try (see Appendix [CJ is still 
not sufficiently good to qualitatively reproduce the exact 
spectral functions in the limit of very strong electron- 
plasmon coupling. 

The fully self-consistent calculations with GT^W P 
seem to only marginally affect the lineshape of the plas- 
mon side-band peaks in comparison to their non self- 
consistent counterpart. 

Note that a fine analysis of the comparison between the 
exact results and the diagrammatic perturbation results 
with vertex correction is difficult to perform in Figs. [5] 
and [3j as the calculations were done for different num- 
bers of w-grid points . It was necessary to perform the 
calculations in that way because the vertex corrections 
scale as N% as shown in Ref. [59j . Therefore we have 
performed the corresponding calculations with a lower 
number of points N u — 1579 for the bottom panels of 
Figs. [5] and [21 instead of N u = 16385 points for the top 
panels, in order to have tractable computational costs. 
Our NEGF code works with a finite broadening related 
to the number of grid points to deal with sharply peaked 
and/or discontinuous functions, hence the different line- 
shape in the spectral functions in the top and bottom 
panels of Figs. [2] and [3] respectively. This numerical ex- 
tra broadening affects only the width of the peaks and 



the global amplitude of the spectral functions, though 
all spectral functions are always normalized. There is no 
major problem with the spectral information contained 
in A(lu). We have discussed in detail the effects of this 
extra broadening in Ref. |59j . 

In addition, we want to add that our results con- 
firm those obtained in earlier studies, see for example 
Refs. [38.39.45.62]. However our self-consistent scheme 
for calculating the second-order diagrams by starting 
with the GW-like Green's function allows us to avoid the 
problem of negative spectral densities (at least within 
the range of parameters we have explored) that were ob- 
tained in Refs. HIE]. 



3. Finite temperatures 

For finite temperatures, the exact result provided by 
Eq. (fT!?)) can be generalized from a thermodynamical 
average over the boson statistics within a canonical 
ensembl o 58 i 61 . In addition to the peaks at e c — nu) p 
(n > 0), one also sees spectral information at e c + nuj p 
(n > 1) which corresponds to absorption of the thermally 
populated plasmons, as shown in Figure 2] 

The results for the spectral functions obtained from the 
diagrammatic expansion of the self-energy and of the ver- 
tex functions as shown in Fig. [1] are shown in Figs. 2] and 
[5] Qualitatively we obtain similar effects of the second- 
order diagrams on the spectral functions as in the case 
of zero temperature. Note that however, for finite tem- 
peratures, the dependence of the lineshape upon the ex- 
tra broadening related to the number of w-grid points 
is much less important, since the thermal broadening 
is dominating. In Fig. 2] we see that, as for the zero- 
temperature case, the self-consistent GW P calculations 
generate the correct series of peaks with the plasmon 
emission sideband peaks again appearing at too low en- 
ergies. However the new plasmon absorption peak just 
above the main peak is almost at the correct energy po- 
sition. 

We do not yet have an accurate explanation for the 
tiny shoulder-like feature around the Fermi level in the 
top panel of Fig. [4] However, this feature is related to 
plasmon absorption processes since at the chosen tem- 
perature the plasmon mode can be thermally populated. 
Nonetheless, it is clear that the feature disappears when 
performing the calculations with an extra broadening 
(i.e. introducing an effective finite lifetime for the plas- 
mon mode). 

When we consider the strong coupling case, shown in 
Fig. [SJ we find that for all levels of approximation the 
lineshape is strongly broadened, washing out most of the 
features. 

We can conclude that, within the limit of the S-model 
and for both the zero-temperature and finite-temperature 
cases, the various GW approximations are much more 
accurate for the high-density regime. For the low-density 
electron gas both the GW peak positions and lineshapes 
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FIG. 4: (Color online) Finite-temperature equilibrium spec- 
tral functions A(u>) for the high-density electron gas with 
rs = 2 and a finite temperature kT — 0.2 corresponding 
to ujp/kT = 3.062. Top panel: Exact results and calcu- 
lations for different self-energies S = GW P and GWjT . 
Bottom panel: Results for different self-energies E = GW P , 
G(F (0) + rf 1 f)W p and G(T (0) + rf^Wp (see Fig. [Q with 
fewer grid points (JV W = 1579), giving an extra broadening. 



are poor in comparison to the exact results, although 
the separation between the plasmon sideband peaks is 
correctly reproduced. 



E. Spectral function of pure jellium and vertex 
corrections 

In this section, we compare different approximations 
for the vertex corrections for another model system: the 
pure jellium model (without a distinct core level). The 
spectral functions in this system are evaluated in the 
zero-temperature limit within conventional Green's func- 



FIG. 5: (Color online) Finite-temperature equilibrium spec- 
tral functions A(u) for the low-density electron gas with 
rs — 4 and a finite temperature kT — 0.2 (corresponding 
to u)p/kT = f.083) and with fewer grid points N u = 1579. 
Calculations for different self-energies E = GW P , GWp w , 
G{T {0) +T^)W P and G{T {0) +Y\^)W p (see Fig.© are shown. 



tions calculations^. 

It is expected from the original work of Hedin et 
al£ and also of Shirley^ 5 - that the exact spectral func- 
tion of pure jellium should show several plasmon res- 
onances below the main quasiparticle peak. However, 
we do not observe any such peaks (see Fig. [5]) when it- 
erating the Green's function to self-consistency within 
the GW approximation, nor when we use model ver- 
tex correction al 45 . These vertex corrections were how- 
ever supposed to provide an exact description of screened 
Coulomb interaction W for the jellium model. 

Any self-consistent iteration has the effect of broaden- 
ing the occupied bandwidth (a feature which is known 
to be unphysical) as evidenced by the shift in the main 
quasiparticle peak at the bottom of the band seen in 
Fig. [6] The model vertex corrections tested do not rem- 
edy this behavior, nor do they lead to any multi-plasmon 
resonances. We consider two different models for the 
vertex corrections: Firstly, a strictly local vertex cor- 
rection applied in the screening, annotated Wo and mod- 
elled directly by the LDA exchange-correlation kernel as 
described by Del Sole et al&. Secondly, the other ver- 
tex correction incorporates a momentum-dependent lo- 
cal field factor modelled on exact quantum Monte Carlo 
results for jellium, as described by Shirley^ 5 - (annotated 
Ws). 

In general, the difference between the two different 
types (static vs. g-dependent) of vertex corrections im- 
plemented is practically negligible in the spectral func- 
tions. This shows that the screened interaction can be 
very insensitive to the exact type of vertex correction 
used, in contrast to the self-energy. With a self-consistent 
calculation, we also observe the broadening of spectral 
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FIG. 6: (Color online) Spectral function A^{uj) at k = and 
for (a) rs = 2 and (b) rs = 4 for the pure jellium model. 
Each set of curves shows a one-shot calculation for Go (black 
data) and self-consistent iterations for G (red data). The 
bottom panel of plots (a) and (b) show standard GoWo and 
GWo, the top panels show GoWo and GWq with a local vertex 
correction in Wo as described in Ref. [f|. GoWs and GWs refer 
to a momentum-dependent vertex correction in Ws defined in 
Ref. [45[ to approximate the exact W of jellium. The results 
including the different vertex corrections are indicated with a 
line (for Wo) and a symbol (for Ws). All chemical potentials 
are aligned at the Fermi energy of the noninteracting gas. 
The positive (negative) deviation of the main quasiparticle 
peak from the origin indicates a narrowing (broadening) of the 
occupied bandwidth. None of these approximations provides 
more than one plasmon satellite in contrast to the expected 
exact result. 



peaks previously noted in Refs. [44ll47l ] . 

This also indicates that the explicit evaluation of the 
second-order diagrammatic vertex correction, ^n), is im- 
perative in order to capture the higher-order plasma sat- 
tclitcs in a metallic system, and in corresponding models 
with a coupling to a core state as shown in the previous 



section. This finding is fully consistent with the previous 
work of Shirley 4 ^ where the vertex function Ttu was ap- 
proximately evaluated within the zero-temperature for- 
malism. 



V. CONCLUSIONS 

We have formally expressed the Hedin's GW equa- 
tions on the Keldysh time-loop contour. This implies 
that within our formalism one can now deal with full 
non-equilibrium conditions for fully interacting electron 
systems. The equilibrium properties of the system are 
obtainable from our formalism as a special case of the 
more general non-equilibrium conditions. 

We have considered in particular the lowest-order ex- 
pansions of the electron self-energy £ and of the vertex 
function T, and compare our results with previous work. 
We have then used our formalism to study a simple model 
of an electron core level coupled to a plasmon mode for 
which exact results for the spectral function are available 
(i.e. the S- model). We have compared our lowest-order 
expansions of the electron self-energy and of the vertex 
function with the exact results, considering the second- 
order diagrams in terms of the plasmon propagator W p . 

We have shown that self-consistent GW-based approx- 
imations (with or without vertex corrections) provide 
a good approximation to the exact results in the limit 
of weak to medium electron-plasmon coupling (i.e. high 
electron-density limit) both at zero and finite tempera- 
tures. Non self-consistent GqW p calculations do not re- 
produce the complete series of plasmon sattelites. How- 
ever the GW based approximations perform quite poorly 
in the strong-coupling limit (i.e. low electron-density 
limit). Vertex corrections generally re- adjust the peak 
positions (the relaxation energy responsible for the renor- 
malization of the core level as well as the plasmon side- 
band peaks) towards the correct result. 

Furthermore we have also analyzed the spectral func- 
tions obtained from conventional equilibrium GW cal- 
culations for the pure jellium model and using different 
approximation for the vertex corrections in W. The cor- 
responding results confirm that the explicit second order 
diagrams for the vertex corrections are needed to obtain 
the full series of plasmon side-band resonances. 

In appendix [Dl we have also addressed an important 
issue about the Dyson-like equation for the time-ordered 
Green's function in the energy represention. We have 
shown that there is a difference between Dyson equation 
for the Green's function obtained at zero-temperature 
and at finite temperature, as already pointed out in 
Ref. Q ■ We have shown that at finite temperature there 
are extra terms in the Dyson equation of the time-ordered 
Green's function. These terms are obtained rigorously 
from the Keldysh time-loop formalism we derived at equi- 
librium, while they were introduced ad hoc by Hedin and 
Lundqvist£ to recover an exact result. 

Finally, we have studied in this paper models of inter- 
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acting electron systems, but we believe that our theoreti- 
cal approach is well-suited for applications towards more 
realistic physical systems, such as the one-dimensional 
plasmon modes recently observed in an atomic-scale 
metal wire deposited on a surface^. 



Appendix C: Lowest order expansion of the vertex 
function T(12; 3) 

1. The rVo) level of approximation: no vertex 
corrections 
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In this section, we derive from our general results the 
more conventional GW approach used in previous stud- 
ies on the ground state properties of molecules, semi- 
conductors, or on the linear response or the full non- 
equilibrium transport properties of nanoscale systems 
driven by an applied external voltag o 28 ' 33 i 34 ' 71 ~ — . 

With no vertex corrections, T(12; 3) is simply given by 
r (0 )(12;3) = <5(12)<5(13). Hence the polarizability P(12) 
and the electron self-energy £(12) are 



or T^^{12). 






X r = x ++ _ x +- 


- = x-+ 


-x~ 


X a = x ++ _ x — 


h = x+~ 


-x~ 


x++ + x— 


- = x+- 


+ x- 


x-+-x+- 


- = x r - 


X a 



Appendix A: Relationship between the different 
Green's functions and self-energies 

The relations between the different components of 
the Green's functions and self-energies on the Keldysh 
time-loop contour are given as usual, with X VlV2 (12) = 



(Al) 



The usual lesser and greater projections are defined 
respectively as X< = X^ and X > = X h , and the 
usual time-ordered (anti-time-ordered) as X t — X ++ 
{X l = X—). 



Appendix B: Rules for analytical continuation 

For the following products Pu-\(t,t') on the time-loop 
contour Ck, 



P(2) = / AB 



C K 



P 



(3) 



ABC 



C K 



P(„) = / A 1 A 2 ...A n , 

Jc K 



we have the following rules for the different components 
Pf^it, t') on the real-time axis: (x = r, a, >, <) 



P |) = J A r B^+A^B a 



P ( < = J A < B a C a + A r B < C a + A r B r C < 

pr _ / AT AT AT pa A a A a A a 

~ J n ' ~ J 



P(12) = -iG(12) G(21), 
S(12) = iG(12) W{21). 



(CI) 



The different components of the polarizability are then 



P§(12) = -iG^(12) G^(21). 



(C2) 



Using Eqs. (|Aip . we hnd that the retarded polarizability 
is given by 



P r (12) = -i[G(12) G(21)] r 

= -iG r (12) G<(21) - iG<(12) G Q (21), 

and the electron self-energy by 



(C3) 



S<(12) = iG<(12) W > (21), 

E r (12) = i[G(12) W{2l)] r (C4) 
= iG r (12) TU < (21) + iG < (12) W a {2\). 

Using the symmetry relations for W and Eqs. (jAip . we 
can easily recast the above equations in the following 
form 



£<(12) = iG<(12) 1U<(12) 

S r (12) = iG''(12) W > {12) + iG<(12) TU r (12). 



(C5) 



These expressions for S and P are just the equivalent of 
Eqs. (3-8) in Ref. [33[ and are similar to the correspond- 
ing expressions in Refs. miHEi!] ' 



2. The Tj!) level of approximation 

With the series expansion T(12;3) = J2 n (12; 3), 
in which the index n represents the number of times the 
screened Coulomb interaction W appears explicitly in the 
series, we take for F(!)(12; 3) 



r (1) (12;3) = j d(4567) 



<5S(12) 
<5G(45) 



G(46)G(75)r(67;3), 



(C6) 
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wherer(67;3) = r (0) (67;3) = 5(67)5(63) and E = iGW. 
Hence r (1) (12; 3) = iW(21) G(13) G(32). 

In the following, we derive the part of the electron self- 
energy and the part of the polarizability arising from Tm 

only. In principle, the full E and P should be calculated 
by using T = + T^). We find for the electron self- 
energy (defined on the contour Ck )'■ 

£(12) = i j d(34) G(13) r (1) (32;4) W(4, 1) 



= i x ij d(34) G(13) W(23) G(34) G(42) W(41). 

(C7) 

The different components E' 71 ' 72 of the self-energy on the 
time-loop contour (with 7/1.2 = ±) are then given by 

E" 1 711 (12) = - ^ 773774 / d(34) G mm (13) W^ 2 " 3 (23) 
mm 

C 3I ' 4 (34) G' 74r ' 2 (42) W mm (4l). 

(C8) 

This self-energy corresponds to the so-called double- 
exchange diagram. Note that we have studied the effects 
of such a diagram in the different context of a propagat- 
ing electron coupled to a local vibration mode, in which 
the bosonic propagator W is replaced by a phonon prop- 
agator D^. 

At the level of approximation, we find that the 
polarizability is given by 

P(12) = -i [ d(34) G(13) G(41) r (1) (34;2) 



= J d(34) G(13) G(41) W{AZ) G(24) G(32), 



(C9) 



with components on Ck given by 

pmm (12) = ^ 7/37/4 / d(34)G I ' 1 " 3 (13) G" 4I ' 1 (41) 

mm 

W V4,13 (A3) C 2, ' 4 (24) G mV2 {32). 

(CIO) 

Here again, and as well as for the self-energy, the re- 
tarded (advanced) part P' r (12) is obtained from P r = 
P ++ — P H . One can then express P r and P H in a 
more compact form involving only terms like X r ' a '^ (with 
X = G,W). 



Appendix D: Time-ordered Green's functions at 
equilibrium 

In this section we discuss in detail the relation between 
time-ordered Green's function (in energy representation) 
for two temperature limits. Differences are expected to 



arise as shown in Chapter IV. 17. of Ref. @|. We use the 
conventional equilibrium many-body perturbation theory 
(MBPT) to determine the time-ordered Green's function 
G*, and the generalization of the Green's function onto 
the Keldysh time-loop contour at equilibrium to deter- 
mine the counterpart of the time-ordered Green's func- 
tion G++. 

From MBPT, the time-ordered Green's function satis- 
fies the Dyson- like equation G* = g t +g t H t G t and the cor- 
responding time-ordered Green's function obtained from 
the Keldysh time-loop expansion satisfies the correspond- 
ing Dyson-like equation G ++ = g ++ + ((/EG) ++ . In prin- 
ciple, from the conventional definition we have g l = g ++ 
and and should have G* = G ++ . 

It is easy to show that from the rules of analytical con- 



tinuation G + 

C++ = g^ 



; + + . 



(</£G) ++ is expanded as follows 



~G H 



- ,g++E < G > + «?<E > G ++ + ff <E~ G++ 



(Dl) 



and after further manipulation (using the notation 
(g/GY = (g/G)++), 



G t = 5 t +( ff t S t - 5 <E>)G t -( 5 S)<G ; 



(D2) 



So, strictly speaking, the non-equilibrium formalism in- 
troduces two extra terms <? < E > G t and (gE) < G > in the 
Dyson equation for G*. 

We now analyze these two terms in more detail. First 
of all, we recall that at equilibrium or in a steady state, 
the Green's functions and self-energies depend only on 
the time difference of their argument and can be Fourier 
transformed with a single energy argument. We then 
have the following expression 

G» = .g 4 M + (s'MEV) - , 9 <HE>H) G*(u,) 
-( 5 £)<MG>M. 

(D3) 

Furthermore, at equilibrium or in a steady state, the 
lesser and greater components of either a Green's func- 
tion or a self-energy (X>) can be expressed in terms of 
the corresponding advanced and retarded quantity and a 
distribution function^—, i.e. 



(D4) 



At equilibrium f^(u>) = /o*(w) and for a system of fer- 
moins, Jq is given by the Fermi-Dirac distribution func- 
tion / oq (w) = 1/(1 + exp0(uj - /j cq )) and /> = / cq - 1 
(with = l/kT). 

At zero temperature, the Fermi-Dirac distribution 
takes only two different values, / cq = 1 or 0. Hence 
we have the property (/ cq ) 2 = / oq , which implies that 
f<{cj)f>(uj) = / cq (/ cq - 1) = 0. Consequently any 
products of the kind X< {uj)Y> (uj) or X > {u))Y < {u) van- 
ish. Therefore we recover from the Keldysh time-loop for- 
malism Eq. (|D3|) at zero temperature, the conventional 
Dyson equation G* = <?* + g l E* G* as expected. 
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At finite temperature /<(a;)/ > (w) = / cq (/ cq - 1) = 
kTduf ccl 7^ 0, and the product /q/q gives a sharply 
peaked function at the Fermi level fi cq = £f with a width 
of approximately kT. 

We now check the individual contribution of each term 
# < E > and (g'E) < G > , first for a specific case (i.e. the 
quasi-particle approximation) and then for the general 
case. 

In a quasi-particle scheme, i.e. when a single index k is 
good enough to represent the quantum states (with en- 
ergy £k), the Green's functions and the self-energies in 
the absence and in the presence of interaction are diago- 
nal in this representation. We have 



function g l : 

G* = g 1 (1 + (E ( - {g t )- l g<Y>) G l - ( 3 t )- 1 ( 5 E)<G>) , 

(D8) 

with (#E)< = g<Y, a + .g r E<. 

By using the equation of motion of the non-interacting 
time-ordered Green's function g t : 



ij--/ l0 (l)) 9 i (12) = 5(12), 



(D9) 



it is straightforward to find that 
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(<? t )- 1 (13)= (i^-- VI)) 5(13). (DID) 



-foiM - flfc)(w) x -fo,k&t - E °)M (D5) and consequently 



For purely fermionic systems at equilibrium, one 
usually has 3mE^(/i cq ) = Oii, and therefore 
9k i^ eq )^k (A* eq ) — 0. When SSmEjr also vanishes in 
the energy window around the Fermi level, defined by 
fafo 7^ 0, then the product g k (ui)T, k (lo) also vanishes. 
When there are no eigenvalues Ek (of the non-interacting 
system) within this energy window, then once more we 
have ff <H£>H~0. 

Otherwise g k (w)E^ - (cj) = Z k 5{u) — Ek) with Zk = 
4 7 r(/ </ >3mE£M) a)=£fc; 

For the second correction term, we have 

{g^)k{u)<G>{uj) 

= -tfMgWu ~ m%)(<») x -foAGl - GS)(«) 



For the quasi-particle scheme, 3mE^° ~ ±vq around 
the Fermi level /x eq ± kT, and we find that 

( 5fc E fc )<G> = -4Trf<f>Z k $leY. r k (e k ) 5(lo - e k )S(cu - e k ) 

(D7) 

with Z k 1 = 1 — (d'SieT, r k / du})uj = i k being the effective mass 
renormalisation parameter and Ek = Ek + ^teYT k being the 
renormalized eigenvalue. Hence the product (pfcSfe) < G^ 
vanishes because in general one has Ek ^ £k- In the 
opposite case when e k = Ek for some quantum states, the 
product (gk^k) < G k also vanishes because then 5ieE£ = 
0. 

Therefore our analysis show that, in the quasi-particle 
scheme at finite temperature, Eq. (|D3|) reduces to the 
conventional Dyson equation G\ — g\ + g\ E| G\ as 
expected. 

Now we need to check what is happening to the two 
contributions g < E > and (gH) < G > beyond the quasi- 
particle approximation. For that we can proceed fur- 
ther: going back to the full time-dependence of Eq. (|D2|) 
and factorizing the non-interacting time-ordered Green's 



(D6) 



{g t )-'g<{lA)= d3 («;*)-' (1%< (34) 
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(Dll) 



i^--h Q (l))g<(U) = 0, 



the last equality comes from the definition of 
5 ,< (14). Similarly one can find that (g t )~ 1 g r = 
Jd3 G?*)- 1 (13)^(34) = tf(14). 

Hence Eq. (|D8[) is transformed into 



G = g 
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E*G* 



E<G ; 



(D12) 



where the last term E < G > satisfies the detailed balance 
equation at equilibrium-^: E < G > = E > G < . 

Eq. (|D12|) is the most general expression for G* and is 
the most important result of this section. It is interest- 
ing to note that Eq. (|D12[) is the equivalent of Eq. (17.9) 
derived in Ref. [jj. However in our approach, the extra 
term E < G > is obtained rigorously from the use of the 
general Kcldysh time-loop contour formalism. While in 
Ref. 0], Hedin and Lundqvist introduced this correction 
term ad hoc in the Dyson equation for the finite temper- 
ature time-ordered Green's function in order to recover 
the proper limit of the independent particle case. 

Once more one can show that, after Fourier transform- 
ing, the product E < G > vanishes at equilibrium and at 
zero temperature because of Eq. (|D3|1 and f < f > = 0. 
Within the quasi-particle scheme at finite temperature, 
we have £<(w)G>(w) = -4/< fc /> fc 3m££(u,) 3mG£(a;). 
Thus, one needs to check the contributions of the spectral 
information in 3mSj(u) and in 3mGJ(w) (in the energy 
window defined by f < k fQ k around the Fermi level) to see 
if the product E^ G^ vanishes (as shown above). 

We conclude this appendix by saying that there is 
indeed a difference between the Dyson equations for 
the time-ordered Green's functions at zero and finite 
temperatu re 3 ' 69 > 76 . This result by no means contradicts 
the fact that the Green's functions on the Keldysh con- 
tour, the time-ordered Green's function at zero temper- 
ature and the Matsubara temperature Green's function 
of imaginary argument all obey the same formal Dyson 
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equation. Our derivations provide a rigorous mathemati- 
cal result for the finite temperature time-ordered Green's 
function (in the energy representation) which satisfies a 
Dyson equation with an extra term as introduced in an 
ad-hoc way in Chap IV. 17. of Ref. Q. 

In our calculations, the correction term £ < G > is au- 
tomatically taken into account since we work with the 
Keldysh time-loop formalism. We have checked nu- 
merically that the £ < G > indeed vanishes at zero tem- 
perature. For finite temperatures we have found that 
£ < G > ~ in the energy window defined by J^/q ^ 
since most of the spectral weight is far below the Fermi 
level (see Figures [5] to [5]) . However, in the limit of very 
high temperatures (i.e. oj p /kT <C 1), the energy window 



defined by Jq/^ ^ is wide and the product E < G > 
does not vanish; though the corrections are two orders 
of magnitude smaller than the amplitude of the Green's 
function G 4 itself. 

It would be interesting to find real cases of interact- 
ing electron systems (probably of low dimensionality) for 
which the correction term £ < G > is not negligible. At 
finite but low temperatures, systems with a strong spec- 
tral density around the Fermi level (i.e. presenting the 
Kondo effect) at low temperature should be a good ex- 
ample. The high temperature limit for metallic systems 
represents another interesting case as shown, for exam- 
ple, in Ref. [10]. 
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